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INTRODUCTION 


Two  axisymmetric  ground  motion  calculations  were  carried 
out  by  Applied  Theory,  Incorporated  (ATI)  to  simulate  decoupled 
explosions  in  mined  cavities  in  salt.  Results  from  these  cal- 
culations were  provided  to  Systems,  Science  and  Software  (SJ) 
for  analysis  to  determine  the  character  of  the  teleseismic 
signatures  of  these  decoupled  explosions. 

The  ATI  calculations  were  in  a homogenous  whole  space 
and  were  carried  into  the  elastic  regime  where  the  divergence 
and  curl  of  the  displacement  field  were  monitored.  From  these 
quantities  an  equivalent  elastic  source  representation  was  com- 
puted for  each  of  the  two  calculations.  In  the  first  case  the 
mined  cavity  was  spherical  and  the  equivalent  elastic  source  is 
a center  of  dilatation.  In  the  second  calculation  the  cavity 
was  an  ellipsoid  of  revolution.  In  this  case,  the  equivalent 
elastic  source  includes  higher  order  terms  representing  the 
departure  from  spherical  symmetry. 

In  Section  2 we  briefly  describe  the  material  properties 
and  geometry  for  the  two  source  calculations.  In  Section  3,  the 
method  for  determining  the  equivalent  elastic  source  is 
described.  It  is  based  on  an  expansion  of  the  displacement 
field  in  spherical  harmonics.  The  source  is  then  described 
by  a series  of  multipole  coefficients.  An  advantage  cf  this 
representation  is  that  the  near-  and  far-field  terms  and  the 
P and  S waves  are  conveniently  separated.  Also,  the  long- 
period  asymptotic  behavior  of  the  coefficients  can  be  bounded. 

In  Sections  4 and  5,  we  describe  the  multipole  coeffi- 
cients in  the  time  and  frequency  domains.  While  the  time 
domain  coefficients  are  exact,  apart  from  numerical  error, 
there  is  some  ambiguity  in  the  Fourier  transformed  coeffi- 
cients at  long  periods.  This  arises  from  the  fact  that  the 
calculations  are  carried  out  for  only  a few  tenths  of  a second 
while  our  primary  interest  is  in  the  character  of  the  solution 
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*t  periods  of  1 to  30  seconds.  The  ambiguity  results  from  the 
fact  that  we  cannot  be  sure  the  values  at  the  last  computed 
time  step  are  close  to  the  true  static  solution.  This  is 
particularly  a problem  for  the  S wave  portion  of  the  solution. 

The  best  we  can  do  is  to  attempt  to  bound  the  contribution  of 
the  long  period  S waves. 

As  far  as  the  spherically  symmetric  portion  of  the  field 
is  concerned,  we  find  that  the  spherical  cavity  is  slightly  more 
decoupled  (20  to  25  percent)  than  the  elliptical  cavity.  The 
higher  order  terms  increase  the  radiation  from  the  latter,  in- 
creasing the  decoupling  difference  between  the  two.  Comparing 
other  estimates  for  the  source  function  for  a tamped  explosion 
in  salt  we  find  decoupling  factors  ranging  from  300  at  long  period 
to  less  than  100  at  frequencies  greater  than  3 Hertz. 

Using  the  equivalent  elastic  source  representation, 
body  and  surface  waves  synthetic  seismograms  were  computed  for 
the  two  sources.  Results  are  presented  in  Sections  6 and  7. 

While  changing  interference  between  the  direct  P,  pP  and  sP 
phases  can  lead  to  minor  variations  in  m^  and  the  waveforms, 
we  conclude  that  radiation  pattern  affects  are  unimportant 
for  short  period  body  waves  and  m^.  That  is,  the  m^  for  the 
ellipsoidal  cavity  is  not  substantially  different  than  that 
for  the  spherically  symmetric  cavity.  For  surface  waves  the 
results  are  very  much  the  same.  It  is  true  that  the  ellipsoidal 
cavity  will  generate  Love  waves.  However,  the  Rayleigh  wave 
radiation  patterns  do  not  differ  substantially  from  the  spheri- 
cally symmetric  case.  These  conclusions  are  based  on  a series 
of  calculations  assuming  the  largest  reasonable  contribution 
from  the  S wave  portion  of  the  field  and  rotating  the  source 
to  orientations  expected  to  give  optimum  radiation  pattern 
effects. 


In  Section  8 we  compare  the  results  obtained  here  to 
results  from  somewhat  similar  calculations  done  several  years 
ago  at  S*  and  reported  by  Cherry,  Bache  and  Patch  (1975)  . 
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These  calculations  were  for  tamped  explosions  in  cylindrical 
cavities  and  lead  to  sources  with  radiation  patterns  similar 
to  those  from  the  ellipsoidal  cavity.  For  these  earlier 
calculations  the  results  were  that  the  m^  and  Ms  were  not 
substantially  different  than  for  the  same  sized  explosion 
in  a spherical  cavity. 


SOURCE  CALCULATIONS 


Two  axisymmetric  ground  motion  calculations  were  carried 
out  by  ATI  to  simulate  explosions  in  mined  cavities  in  salt. 

One  cavity  was  spherical  and  the  other  was  an  ellipsoid  of 
revolution  or  prolate  spheriod.  Pertinent  information  provided 
to  S1  about  these  calculations  is  summarized  in  this  section. 

The  two  cavities  were  assumed  to  be  mined  from  a medium  in  a 
uniform  hydrostatic  stress  of  181  bars.  The  calculations  were 
initiated  with  an  initial  stress  field  that  included  the  effect 
of  mining  the  cavities.  The  spherical  cavity  has  a radius  of 
66  meters  while  the  prolate  spheriod  has  a 3/1  aspect  ratio. 

The  volume  of  the  two  is  the  same.  The  pertinent  dimensions 
and  material  properties  are  listed  in  Table  1. 

The  calculations  were  carried  into  the  elastic  regime 
and  the  divergence  and  curl  of  the  displacement  field  were 
monitored  at  a constant  radius  of  928  meters.  These  data,  the 
time  histories  of  the  divergence  and  curl  at  one  degree  intervals 
on  the  radius,  were  provided  to  Ss  for  analysis.  The  remainder 
of  this  report  is  devoted  to  describing  the  results  of  that 
analysis. 


EQUIVALENT  ELASTIC  SOURCE  - THE  METHOD 


This  problem  is  quite  similar  to  those  studied  by  Cherry 
Bache  and  Patch  (1975)  who  were  interested  in  tamped  explosions 
in  cylindrical  cavities.  The  methods  applied  are  identical  to 
those  in  the  earlier  work.  The  geometry  is  shown  in  Figure  1. 
We  first  compute  multipole  coefficients  in  the  time  domain  from 


TABLE  1 


PROPERTIES  FOR  ATI  SOURCE  CALCULATIONS  IN  SALT 


Depth  of  Burial  (m) 

819 

Density,  p (gm/cm3) 

2.24 

P Wave  Velocity,  a (km/sec) 

4.4152 

S Wave  Velocity,  S (km/sec) 

2.5 

Yield  (KT) 

25 

Initial  Cavity  Volume  (m3) 

1.2  x 10 

Spherical  Cavity  Diameter  (m) 

132 

Prolate  Spheroid  Major  Axis  (m) 

275 

Prolate  Spheroid  Minor  Axis  (m) 

92 
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where  the  curl  and  divergence  are  monitored  on  the  spherical  sur- 
face with  radius  R.  The  axisymmetry  is  accounted  for  in  deriving 
Equation  (1). 


The  next  step  is  to  compute  the  frequency  domain  multipole 
coefficients  from  the  following  transformations 
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ft 


a - 1,  2,  3,  4,  with  a - 4 indicating  the  dilatational  or  P wave 
coefficients.  Having  these  coefficients,  we  can  compute  the  whole 
space  displacements  from 


u(R,w) 


-T-  ?x4  + -K  v X X 

V V * 


(5) 


where  k and  kfl  are  the  P and  S wave  numbers. 

a d 

Using  the  multipole  coefficients  to  represent  the  source, 
we  can  compute  body  and  surface  waves  in  layered  earth  models 
using  methods  outlined  by  Bache  and  Harkrider  (1976)  and  Harkrider 
and  Archambeau  (1978).  While  the  original  coefficients  are  com- 
puted with  the  source  oriented  as  shown  in  Figure  1,  they  can  be 
rotated  (Minster,  1976)  to  represent  a source  at  an  arbitrary 
orientation  with  respect  to  a fixed  coordinate  system. 

The  spherically  symmetric  calculation  can  also  be  repre- 
sented by  a reduced  displacement  potential  defined  by 


u(R, t) 


3 

JR 


y (t-R/q) 
R 


(6) 


In  terms  of  the  multipolar  expansion  we  have 


»£'  <•> 


I .2  A 

~ — k y (u >)  , 

a 


(7) 


and  all  other  multipoles  vanish.  Also, 


£ , .,2  — ik  R 

y(w)  - --tt--  e a - 


b) 


X4(R»w) 


(8) 


where  x4  is  simply  the  Fourier  transform  of  the  divergence. 
Another  useful  relationship  is 


Ra 


- 7 • u(R, t) 
% 


(9) 
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From  these  equations  we  can  derive  some  limiting  behavior 
for  the  multipole  coefficients  that  can  be  used  to  check  the 
solution.  We  know 

❖ 

lim  (uj)  - lim  ’i'(t)  - C (10) 

oj-*o  t-*-® 


where  C is  some  constant.  Then  from  Equations  (9)  and  (10) 


lim  / 7 • u(R,t)  dt  ■ 0 

■ o 


(11) 


From  Equation  (1)  we  see  that 

(R,t)  - - 7 • u(R,t)  (12) 

and  the  limiting  relationship  in  Equation  (11)  applies  to  the 
time  domain  monopole  as  well. 

For  the  higher  order  coefficients  we  have  only  the 
requirement  that 

lim  (R,t)  - C,  n > 2 , (13) 

t-*®  “II 


and  C may  be  a positive  or  negative  constant  or  zero.  If 
C j*  0,  then 

Um  A^“)  (io)  a u>n  , (14) 


otherwise  the  falloff  is  at 


n+1 


CO 


or  steeper  at  low  frequency. 


9 


4 


MULTIPOLE  COEFFICIENTS  IN  THE  TIME  DOMAIN 


r 
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The  first  step  is  to  apply  Equation  (1)  and  obtain  the 
multipole  coefficients  in  the  time  domain.  The  monopole  and 
quadrupole  (double-couple)  terms  for  the  spherical  cavity  cal- 
culation are  shown  in  Figure  2.  If  the  calculation  were  spher- 
ically symmetric,  as  it  would  be  in  the  absence  of  numerical 
errors  in  the  finite  difference  code,  the  quadrupole  terms 
would  vanish.  Indeed,  we  see  that  the  peak  values  are  down  by 
one  or  two  orders  of  magnitude.  Also,  note  from  Equation  (12) 
that  the  monopole  is  just  the  negative  of  the  divergence  of  the 
displacement  field  for  a spherically  symmetric  explosion.  The 
higher  order  terms  are  nonzero  because  the  V • u is  not  exactly 
the  same  at  all  points  monitored  in  the  ATI  calculations. 

The  time  domain  multipole  coefficients  for  the  ellipsoidal 
cavity  are  shown  in  Figure  3.  The  coefficients  for  n <_  4,  the 
monopole,  quadrupole  and  octupole,  are  plotted.  Comparing  to  the 
coefficients  in  the  spherical  cavity  case,  we  see  that  the 
monopole  is  somewhat  smaller  but  is  otherwise  not  greatly  dif- 
ferent. However,  the  higher  order  terms  are  now  the  same  order 
of  magnitude  as  the  leading  term,  the  monopole.  It  is  difficult 
to  estimate  the  relative  importance  of  these  coefficients  because 
of  differences  in  frequency  content.  For  this  reason  the  radia- 
tion field  will  be  studied  in  the  frequency  domain. 

A prominent  feature  of  the  P wave  coefficients  is  an  oscil- 
lation at  30  Hz.  This  is  apparently  related  to  a natural  frequency 
of  the  numerical  grid  (Neil  Perl,  ATI,  personal  communication) . 

5.  MULTIPOLE  COEFFICIENTS  IN  THE  FREQUENCY  DOMAIN 

To  obtain  the  frequency-domain  multipole  coefficients 
we  apply  the  transformation  in  Equation  (3).  While  the  time 
domain  coefficients  of  Figures  2 and  3 are  exact,  apart  from 
numerical  error,  there  is  some  ambiguity  in  the  transformed 
coefficients  at  long  periods.  Examining  the  coefficients  for 
the  ellipsoidal  cavity,  Figure  3,  the  ambiguity  results  from 
the  fact  that  we  cannot  be  sure  the  values  at  the  last  time 
step  are  close  to  the  static  solution. 


10 


I 


M ) 

The  dilatational  coefficients,  A , have  a signal  lasting 

no 

about  0.38  seconds.  Therefore,  estimates  of  these  coefficients 
at  frequencies  less  than  2.6  Hz  will  be  limited  by  the  accuracy 
of  our  estimate  of  the  static  solution.  For  the  rotational 
coefficients  the  same  is  true  for  frequencies  less  that  4 Hz. 


s 
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(4) 

For  obtaining  (w)  we  take  advantage  of  the  limiting 

relationships  of  Equations  (11)  and  (12) . The  quantity 

g (T)  - / T aIq5  (R,t)  dt 
Jo 

is  plotted  in  Figure  4 for  the  two  source  calculations.  Since 

g(T)  is  nearly  zero  at  long  times,  we  can  with  confidence  set 

g(T)  ■ 0 for  T greater  than  the  computational  time.  Then  our 

(4 ) 

long  period  estimate  of  AqQ' (<d)  should  be  quite  good. 

For  the  quadrupole  and  higher  order  terms  the  long 
period  values  depend  on  our  estimate  of  the  static  values. 
There  are  many  ways  this  estimate  could  be  made,  but  the 
simplest  choices  are  these: 

1.  Set  the  coefficients  to  zero  at  times 
greater  than  the  last  time  step. 

2.  Assume  the  value  at  the  last  time  step 
is  the  static  value. 

As  is  clear  from  the  plots  in  Figure  3,  neither  choice  is 
correct,  especially  for  the  octupole  terms.  However,  as 
was  explained  in  connection  with  Equation  (14) , if  we  choose 

the  first  alternative,  the  coefficients  A^  will  fall  off  as 

n+1  nm  n 

u>  at  low  frequency.  Otherwise,  they  will  fall  off  as  <d  . 

Therefore,  we  will  look  at  the  results  from  both  alternatives 

and  attempt  to  bound  the  true  solution. 

The  multipole  coefficients  for  the  ellipsoidal  cavity 

2 

are  plotted  in  Figure  5 after  dividing  by  u>  • It  can  be  shown 


ft 


TIME  (sec) 


Figure  4.  The  time  integral  of  the  monopole  is  plotted 
for  the  two  source  calculations. 
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Figure  5a.  The  "P  wave"  multipole  coefficients  for  the  ellip- 
soidal cavity  are  plotted  after  normalization  by 
w*.  The  arrows  labeled  "Dilatation"  and  "Rotation" 
indicate  the  lowest  frequencies  sampled  by  the 
finite  difference  calculation;  that  is,  the  inverse 
of  the  signal  duration. 


m/(m) 


from  Equations  (4)  and  (5)  that  the  far-field  displacement  is 
directly  proportional  to  a sum  of  terms  of  the  form  /u  . 
Therefore,  the  curves  in  Figure  5 are  proportional  to  the  far- 
field  displacement  spectrum.  From  Figure  5a  we  draw  the  fol- 
lowing conclusions  for  the  dilatational  or  far-field  P wave 

portion  of  the  field  (A^(u>)). 

no 

• The  dilatational  portion  of  the  radiation 
field  is  nearly  spherically  symmetric  for 
periods  longer  than  5 or  10  seconds. 

• The  quadrupole  term  is  significant  at  fre- 
quencies near  one  Hertz,  the  important 
frequencies  for  teleseismic  short  period 
body  waves. 

• It  does  not  matter  very  much  whether  we  assume 

(4) 

the  late-time  double-couple  (A^g  ) is  zero  or 
equal  to  the  value  at  the  last  computed  time 
step.  This  indicates  that  the  computation  is 
essentially  complete  as  far  as  this  component 
is  concerned. 

• The  octupole  term  is  very  small  if  it  has  a 
zero  static  limit.  If  the  last  time  point 
were  the  static  limit,  the  octupole  term  makes 
a large  contribution  near  one  Hertz.  However, 
looking  at  Figure  3,  this  seems  to  be  an  unrea- 
sonable assumption. 

For  the  rotational  or  far-field  S wave  portion  of  the 
radiation  field,  the  story  is  somewhat  different.  We  note  the 
following: 

(2) 

• The  octupole  term,  A^  , is  not  important. 

• For  the  double-couple  term  our  assumption  about 
the  long  time  behavior  makes  a huge  difference 
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at  long  period.  This  indicates  that  the  calcu- 
lation has  not  been  carried  out  long  enough  for 
the  rotational  part  of  the  field  to  pass  completely 
through  the  monitoring  surface.  This  is  also 
apparent  from  the  plots  in  Figure  3. 

• At  the  M period  (a:  20  seconds)  the  S wave 
s 

coefficient  is  either  0.25  or  6.5  times  the 
P wave  coefficient.  The  range  is  therefore 
a factor  of  26.  The  truth  probably  lies 
somewhere  between. 

In  Figure  6 we  compare  the  monopole  for  the  ellipsoidal 
cavity  to  that  for  the  spherical  cavity.  The  spherically  sym- 
metric contribution  to  the  displacement  field  is  smaller  for  the 
spherical  cavity  below  10  Hertz.  This  result  seems  a bit  sur- 
prising because  the  source  yield  and  cavity  volume  are  the  same 
for  both  calculations  and  the  field  for  the  ellipsoidal  cavity 
is  augmented  by  contributions  from  higher  order  terms . 

Radiation  patterns  showing  the  excitation  of  far-field 
P and  S waves  at  selected  frequencies  are  shown  in  Figure  7 for 
the  spherical  cavity  and  Figure  S for  the  ellipsoidal  cavity. 

The  axis  of  symmetry  for  the  latter  is  the  Z-axis  noted  in  the 
figure.  Thus  there  is  no  azimuthal  (about  Z)  dependence  and 
the  patterns  shown  completely  describe  the  radiation.  The 
patterns  for  the  spherical  cavity  show  the  small  deviation 
from  spherical  symmetry  that  was  previously  pointed  out. 

Radiation  patterns  are  shown  in  Figure  8 for  the  two 
alternate  assumptions  for  the  static  behavior  of  the  double- 
couple terms.  In  either  case  the  P wave  pattern  does  not 
deviate  strongly  from  spherical  symmetry  until  we  reach  fairly 
high  frequencies.  Further,  at  1 and  3 Hz  there  is  not  much 
difference  between  the  two  solutions.  However,  at  the  fre- 
quencies important  for  far-field  surface  waves  and  M , 0.05 
and  0.1  Hz,  the  S waves  range  from  very  large  to  very  small, 
depending  upon  our  assumptions  about  the  static  limit. 
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Figure  7 


Radiation  patterns  are  shown  at  six  frequencies  for 
the  spherical  cavity.  The  displacement  amplitude  for 
the  innermost  circle  is  shown  at  the  right,  normalized 
to  that  for  the  0.05  Hz  case. 


The  source  for  the  spherical  cavity  can  also  be  repre- 
sented in  terms  of  a reduced  velocity  potential  using  Equation 
(7)  or  (9).  In  Figure  9 we  show  the  l^l  computed  from  the  7 • 
monitored  at  two  stations  located  at  0 ■ 10#  and  0 ■ 30*.  The 
differences  between  the  two  are  expected  in  view  of  the  fact 
that  this  source  is  not  quite  spherically  symmetric. 


Also  plotted  in  Figure  9 is  the  |$|  computed  from  the 
monopole  (Figure  6)  using  Equation  (7).  The  potentials  com- 
puted directly  from  the  7 • u compared  to  that  from  the  mono- 
pole  differ  by  a significant  amount  at  the  long  period  end. 
Since  the  monopole  represents  a spatial  average  of  the  7 • u; 
that  is. 


(R,t) 


u (R, t ) sin  0d0 


(15) 


from  Equation  (1) , the  estimate  based  on  the  monopole  is 
probably  more  accurate. 

As  an  estimate  of  the  amount  by  which  the  explosions 
studied  here  are  decoupled  in  Figure  10,  we  compare  the  mono- 
pole  source  function  from  Figure  9 to  two  estimates  for  the 
source  function  for  coupled  explosions.  The  Mueller/Murphy 
source  function  (Mueller  and  Murphy,  1971)  is  based  on  a fit 
to  the  SALMON  data.  Therefore,  it  can  be  taken  to  represent 
our  empirical  experience  in  this  material.  The  S1  source 
function  255  was  computed  with  one-dimensional  finite  dif- 
ference methods  for  a tamped  explosion  at  a depth  of  500 
meters  and  is  described  by  Bache,  Cherry  and  Mason  (1976). 

It  represents  computational  results  for  explosions  in  this 
material . 

The  decoupling  factors  are  indicated  on  Figure  10. 

The  decoupling  is  largest  at  long  periods  and  decreases 
rapidly  with  increasing  frequency  above  one  Hertz. 
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Pigure  9.  Three  estimates  of  the  reduced  velocity  potential  for  the  spherical 
cavity  calculations  are  compared. 
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SOURCE 


SOURCE  255 
(Bache,  Cherry  and 
Mason*  1976) 


FREQUENCY  (Hz) 
Decoupling  Factors 


Mueller/Murphy  Source 
Source 255 


The  source  function  from  Figure  9 is  compared  to 
two  estimates  for  the  source  function  for  a fully 
coupled  explosion  in  salt.  The  decoupling  factor 
is  computed  with  the  Mueller /Murphy  source. 


6.  BODY  WAVE  SEISMOGRAMS 

Using  the  multipole  coefficients  presented  in  the  previous 
section  we  can  compute  body  and  surface  wave  seismograms  in  the 
same  way  as  done  by  Bache,  Masso  and  Mason  (1977).  For  the  body 
waves  we  did  the  following: 

• The  source  is  embedded  in  the  crustal  structure 
tabulated  in  Table  2. 

• The  calculations  are  for  a depth  of  820  meters. 

• The  body  waves  were  computed  for  a ray  param- 
eter of  0.079  sec/km.  This  corresponds  to  a 
distance  of  36°.  The  takeoff  angle  for  the 
source  is  then  20.4°  from  the  vertical  for  the 
P and  pP  phases  and  11.4°  from  the  upward  ver- 
tical for  the  sP  phase.  The  delay  times  are 
P-pP  * 0.45  seconds  and  P-sP  ■ 0.71  seconds. 

-4 

• The  P waves  were  multiplied  by  0.95  x 10  to 
represent  the  geometric  spreading  for  a range 
of  36° , 1.77  to  account  for  the  free  surface 
at  the  receiver  and  an  operator  for  anelastic 
attenuation  with  t*  * 0.7. 

• The  response  of  a short  period  seismograph  was 
also  included. 

• The  is  then  computed  from 


“b  " I©5  t 


+ 3.32 


The  synthetic  seismograms  are  shown  in  Figure  11.  From 
top  to  bottom,  the  sources  are  as  follows: 

1.  The  source  255  of  Bache,  Cherry  and  Mason 
(1977)  is  used. 
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2.  The  source  is  the  monopole  for  the  spherical 
cavity. 

3.  The  source  is  the  V(u»)  computed  from  the  V • u 
for  the  spherical  cavity  at  8 ■ 10°. 

The  source  functions  for  the  first  three  seismograms  were  com- 
pared in  Figures  9 and  10.  The  other  seismograms  are  for  the 
elliptical  cavity  for  which  the  multipole  coefficients  were 
shown  in  Figure  5.  These  are: 

4.  The  source  is  the  monopole  term  of  the 
expansion  of  the  radiation  field;  that 
is,  the  spherically  symmetric  part. 

5.  Terms  up  to  n * 4 are  included  with  a zero 
static  value  assumed  for  the  quadrupole. 

The  1 Hz  radiation  pattern  for  this  source 
is  shown  on  the  left  side  of  Figure  8. 

6.  The  source  is  the  same  as  for  Item  5 but 
is  rotated  to  maximize  the  contribution 
from  the  sP  phase. 

7.  Terms  up  to  n ■ 4 are  included  with  a con- 
stant static  value  assumed  for  the  quadru- 
pole. The  source  is  then  rotated  as  for 
Item  6;  that  is,  to  maximize  the  sP  phase. 

The  1 Hz  radiation  pattern  for  this  source 
is  shown  in  Figure  12.  It  should  be  com- 
pared to  the  patterns  in  Figure  8. 

8.  The  source  is  the  same  as  that  shown  in 
Figure  12  for  Item  7 but  has  been  rotated 
54°  counterclockwise  to  give  an  sP  phase 
that  is  on  the  lobe.  That  is,  the  sP 
phase  for  this  case  is  the  negative  of 
that  for  Seismogram  7. 
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Figure  12.  The  radiation  pattern  is  plotted  for  the  source  for 
seismogram  7 of  Figure  11.  The  axis  of  symmetry 
has  been  rotated  38*  from  the  vertical. 
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The  seismograms  for  the  eight  sources  are  shown  in 
Figure  11  together  with  the  m^  data  for  each  one. ' The  results 
are  consistent  with  our  expectations  from  the  analysis  of  the 
multipole  coefficients  and  source  radiation  patterns.  The 
important  points  are  as  follows: 

• Our  tamped  source  255  leads  to  m^  values 
that  are  more  than  two  units  larger  than 
for  the  decoupled  sources. 

• The  three  calculations  for  spherically 
symmetric  sources,  2,  3,  and  4,  give 
about  the  same  result.  The  spherically 
symmetric  part  of  the  field  from  the 
elliptical  cavity  is  slightly  larger 
than  that  from  the  spherical  cavity. 

• The  assumption  about  the  late-time  behavior 
of  the  quadrupole  has  no  importance  for 
short  period  body  waves. 

• Changing  interference  between  P,  pP  and  sP 
as  the  source  orientation  is  varied  leads 
to  some  minor  variations  in  m^  and  the 
waveforms.  This  interference  also  depends 
on  source  depth,  as  is  also  the  case  for 
spherically  symmetric  explosions. 

We  conclude  that  radiation  pattern  effects  are  unimportant 
for  the  short  period  body  waveforms  and  m^. 


7.  SURFACE  WAVE  SEISMOGRAMS 

The  surface  wave  calculations  are  essentially  the  same 
as  in  Bache,  Rodi  and  Harkrider  (1978)  and  Bache,  Masso  and 
Mason  (1977) . The  path  model  is  the  North  American  crustal 
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model  of  McEvilly  (1964) . In  the  source  region  this  is  modified 
by  replacing  the  top  three  kilometers  with  the  top  three  kilo- 
meters of  the  model  in  Table  2.  The  Q operator  is  that  of 
Tryggvason  (1965).  An  LRSM  long  period  seismometer  response  is 
used.  The  seismograms  were  computed  at  3000  km.  The  was 
computed  using  the  Marshall  and  Basham  (1972)  formula  which  is 

M - log  A + 1.08  + P (T) 

9 

1 

where  A is  the  maximum  peak-to-peak  amplitude  near  20  seconds  and 
P (T)  is  a period  correction  tabulated  by  Marshall  and  Basham 
(1972).  The  P (T)  is  0.05  at  24  seconds  and  zero  at  20  seconds. 

The  surface  wave  amplitudes  will  be  proportional  to  the 
long  period  level  of  the  source  functions  discussed  in  Section 
5.  The  only  aspect  that  is  not  entirely  clear  is  the  extent 
of  the  radiation  pattern  effects  introduced  by  the  higher  order 
terms. 

Surface  wave  seismograms  are  shown  in  Figure  13.  The 
sources  for  these  calculations  are  listed  below.  The  first 
three  are  for  spherically  symmetric  sources  and  the  others  are 
for  the  elliptical  cavity. 

1.  S * calculated  source  255  for  a 25  kt  tamped 
explosion  in  salt  is  used, 

2.  The  source  is  the  monopole  for  the  spherical 
cavity . 

❖ 

3.  The  source  is  the  'P  (u»)  computed  from  the 

V • u at  9 * 10°  for  the  spherical  cavity. 

4.  The  source  is  the  monopole,  the  spherically 
symmetric  term,  in  the  expansion  of  the 
radiation  field  for  the  elliptical  cavity. 


(17) 
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Figure  13.  Theoretic  seismograms  are  compared  for  seven  sources. 
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5.  Terms  up  to  n ■ 4 are  included  with  a zero 
static  value  assumed  for  the  quadrupole. 

The  source  orientation  is  vertical. 

6.  Terms  up  to  n ■ 4 are  included  with  a 
constant  static  value  assumed  for  the 
quadrupole.  The  source  orientation  is 
vertical. 

7.  The  source  is  the  same  as  in  Item  6 but 
has  been  rotated  to  the  orientation  shown 
in  Figure  12. 

The  M values  are  indicated  on  the  seismograms  of  Figure 
s 

13.  They  are  entirely  consistent  with  the  spectral  amplitudes 

and  radiation  patterns  shown  in  Section  5.  For  the  spherical 

sources  the  M values  are  as  expected  from  Figures  9 and  10. 

8 

The  source  255  gives  an  Mg  that  is  2.4  to  2.7  units  larger  than 
for  the  decoupled  explosion  in  the  spherical  cavity.  The  multi- 
pole coefficient  source,  the  monopole,  is  0.28  M units  larger 

s 

than  that  estimated  from  a single  V • u time  history.  Also,  the 
spherically  symmetric  part  of  the  source  for  the  elliptical 
cavity  is  somewhat  bigger,  about  0.11  M units,  then  that  for 

S 

the  spherical  cavity. 

The  Mg  values  for  the  seismograms  of  Figure  13  and  a 
number  of  others  are  summarized  in  Table  4.  Also  shown  in 
the  table  are  the  spectral  amplitudes  at  20  seconds.  For  the 
elliptical  cavity  both  the  M and  the  logarithm  of  the  20 
second  spectral  amplitudes  are  normalized  to  the  value  for 
the  monopole  term  taken  alone.  These  values  demonstrate  clearly 
that  the  addition  of  higher  order  terms  to  the  source  repre- 
sentation has  only  a small  effect. 
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8.  COMPARISON  TO  RESULTS  FROM  SIMILAR  CALCULATIONS 

^ Several  years  ago  S 1 studied  the  teleseismic  field  from 

explosions  in  cylindrical  cavities.  The  results  were  reported 
by  Cherry,  Bache  and  Patch  (1975) . The  explosions  were  tamped, 
so  the  decoupling  issues  raised  here  were  not  considered. 

, However,  the  radiation  field  did  include  a large  S wave  compo- 

nent with  the  same  character  as  that  for  the  elliptical  cavity 
studied  here.  Therefore,  it  is  useful  to  compare  results  of 
the  two  studies.  We  should  also  mention  a unique  aspect  of  the 

* work  described  by  Cherry,  Bache  and  Patch.  That  is,  some  small- 
scale  experiments  were  done  to  help  verify  the  computational 
results. 

Two  calculations  were  done.  In  the  first  calculation 
a step  pressure  was  introduced  into  a two-to-one  aspect  ratio 
cylindrical  cavity  as  shown  in  Figure  14.  The  material  response 
was  restricted  to  the  linear  elastic  regime.  The  resulting  far- 
field  radiation  pattern  is  shown  at  several  source  orientations 
in  Figure  15.  This  pattern  applies  for  all  frequencies  below 
about  3 Hz.  Comparing  to  Figure  8,  we  see  it  is  about  the  same 
as  the  pattern  at  f * 0.1  Hz  for  the  zero  static  quadrupole 
case. 

The  second  calculation  studied  by  Cherry,  Bache  and 
Patch  was  a full-scale  simulation  of  a 27  kt  nuclear  explosion 
in  a 3 meter  diameter,  200  meter  long  tunnel  in  saturated  tuff. 

* The  calculation  was  done  by  linking  an  Eulerian  hydrodynamic 
program  with  a Lagrangian  wave  propagation  program.  Two  far- 
field  radiation  patterns  for  this  source  are  shown  in  Figure  16. 
At  long  periods  the  P wave  field  is  significantly  perturbed  but 

* the  S wave  excitation  is  small. 

In  the  previous  report  we  concluded  that  the  short  period 
waveforms  and  m^  were  essentially  the  same  as  from  a spherical 

0 tamped  explosion.  However,  we  reported  a substantial  effect  on 

the  Rayleigh  wave  amplitude  and  M . This  conclusion  was  based 
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Figure  14. 


The  coordinate  system,  geometry  and  physical 
parameters  characterizing  the  linear  elastic 
source  calculation.  The  inset  shows  the  rela- 
tion between  the  cartesian  (1,2,3),  CRAM  (y,z) 
and  spherical  (r,0,$)  coordinate  systems. 


P (Monopol*) 


Figure  15. 


The  far- field  displacement  radiation  pattern 
(f  < 3 Hz)  10  km  from  the  source  of  Figure  14. 
The  P wave  radiation  from  the  monopole  alone  is 
also  indicated,  as  is  the  ratio  of  S to  P wave 
displacements  (on  the  right  half  of  the  figure) 
for  the  total  source. 
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on  calculations  that  were  later  discovered  to  be  in  error, 
this  report  we  will  give  the  correct  results. 


In 


The  Rayleigh  wave  excitation  was  studied  at  three  source 


orientations*  vertical*  horizontal  and  at  45°.  The  M values 

S 


were  calculated  in  the  usual  way.  The  "Mg  enhancement"  was 


then  defined  as  the  difference  between  the  from  the  total 

s 


source  and  from  the  monopole  term  taken  alone.  For  the  source 


oriented  vertically  (as  in  Figures  15  and  16)  the  Mg  enhance- 


ment is  0.09  for  the  linear  elastic  source  and  -0.10  for  the 
27  kt  explosion  in  tuff. 


The  M enhancement  was  calculated  at  a number  of  azimuths 

9 


to  define  the  radiation  pattern  for  the  source  oriented  hori- 
zontally and  at  45°  from  the  vertical.  These  patterns  are 
plotted  in  Figures  17  and  18.  The  difference  between  the 


largest  and  smallest  Mg  enhancement  is  shown  as  well  as  the 


average  over  the  pattern.  It  is  clear  that  the  effect  is 
quite  small.  We  must  conclude  that  this  kind  of  source  asym- 
metry has  only  a weak  potential  effect  on  the  Rayleigh  wave 
amplitudes. 


(a)  One  second  radiation  pattern 


/m 


(b)  Twenty  second  radiation  pattern. 


Figure  16. 


Far-field  displacement  radiation  pattern  10  km  from 
the  nonlinear  source  of  Section  7.  Both  1 and  20 
second  patterns  are  shown  and  the  ratio  of  S to  P 
wave  displacements  are  indicated  on  the  right  half 
of  each  pattern.  The  P wave  radiation  from  the 
monopole  alone  is  also  indicated. 
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The  azimuthal  distribution  of  Ms  enhancement  is 
plotted  for  the  27  kt  explosion  calculation 
oriented  horizontally  (left)  and  at  45*  (right) 
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